Nuclear structure of ^"S and its implications for nucleosynthesis in classical novae 
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Background: The uncertainty in the ^^P(p, y)"'''S reaction rate over the temperature range of 0.1 < T < 
1.3 GK was previously determined to span ~4 orders of magnitude due to the uncertain location of two previ- 
ously unobserved 3+ and 2+ resonances in the Ey = 4.7 - 4.8 MeV region in ■'"S. Therefore, the abundances 
of silicon isotopes synthesized in novae, which are relevant for the identification of presolar grains of putative 
nova origin, were uncertain by a factor of 3. Purpose: (a) To investigate the level structure of -"'S above the 
proton threshold (4394.9(7) keV) via charged-particle spectroscopy using the ^^S(p,t)^^S reaction and in-beam 
y-ray spectroscopy using the Si( He, ny) S reaction to calculate the P(p,Y) S reaction rate, (b) To ex- 
plore the impact of this rate on the abundances of silicon isotopes synthesized in novae. Method: Differential 
cross sections of the ■'^S(p,/)-''*S reaction were measured at 34.5 MeV. Distorted Wave Born Approximation 
(DWBA) calculations were performed to constrain the spin-parity assignments of the observed levels, including 
the two astrophysically important levels. An energy level scheme was deduced from y-Y coincidence measure- 
ments using the ^^Si("'He, ny)^^S reaction. Spin-parity assignments based on measurements of y-ray angular 
distributions and y-y directional correlation from oriented nuclei were made for most of the observed levels 
of 3°S. Results: The resonance energies corresponding to the states with 4.5 MeV < Z:^- < 6 MeV, including 
the two astrophysically important states predicted previously, are measured with significantly better precision 
than before. The spin-parity assignments of both astrophysically important resonances, as well as the existence 
of a previously observed tentative state at ii.v ~ 5.95 MeV are confirmed. The uncertainty in the rate of the 
^^P(P,Y)^''S reaction is substantially reduced over the temperature range of interest. Finally, the influence of 
this rate on the abundance ratios of silicon isotopes synthesized in novae are obtained via ID hydrodynamic 
nova simulations. Conclusions: The uncertainty in the ^^P(p, y)"'*^S reaction rate is reduced to the point that 
it no longer affects the silicon isotopic abundance ratios significantly, and thus the results of our nova hydro 
simulation for the nucleosynthesis in the Si-Ca mass region are more reliable than before. 

PACS numbers: 26.30.Ca,25.40.Hs,23.20.En,23.20.Lv 



I. ASTROPHYSICAL MOTIVATION 



Classical nova outbursts are caused by explosive hydrogen 
burning as a result of a thermonuclear runaway in the enve- 



*Electronic address: kiana.setoodeh@gmail.com 

^Present Address: Department of Physics & Astronomy, Louisiana State Uni- 
versity, Baton Rouge. LA 70803, USA 

^Present Address: Department of Physics & Astronomy, Michigan State 
University, East Lansing, MI 48824, USA; and National Superconducting 
Cyclotron Laboratory, Michigan State University, East Lansing, MI 48824, 
USA. 



lope accreted from a main sequence star onto a white dwarf 
in a close semi-detached binary system. Simulations [1] show 
that peak temperatures reached in the thermonuclear runaway 
are typically in the 0.1- 0.4 GK range, and the ejecta show 
significant nuclear processing. The dominant nuclear reaction 
flow proceeds close to the valley of stability on the proton-rich 
side and is dominated by a series of (/?, y) and {p, a) reactions, 
as well as p+-decays. Classical nova outbursts are thought to 
be the major source of '^N, '^O and to some extent '^C in the 
Galaxy [2] and contribute to the abundances of other species 
with masses up to A w 40, including ^^Al. 

The ejecta of classical novae are studied by systematic in- 
frared observations [3, 4] which reveal episodes of dust for- 
mation following a nova outburst. Several candidate presolar 
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grains of nova origin have been found [5, 6], most of which are 
of siHcon carbide (SiC) type. These grains show abundance 
anomahes for some isotopes (compared with the average so- 
lar system isotopic abundances), e.g., close to or slightly lower 
than solar ^^Si/^^Si ratios and higher than solar ^"Si/^^^Si ra- 
tios [7]. 

In order to reach a quantitative agreement between the iso- 
topic abundances observed in the presolar grains [6] and those 
predicted by simulations [7], nova nucleosynthesis models re- 
quire some dilution. Thus, the mixing between the material in 
nova ejecta and the solar-like material must be understood to 
tighten the links between nova nucleosynthesis and presolar 
grains. Also, a better knowledge of the rates of the reactions 
that affect nova nucleosynthesis is required to better under- 
stand the origin of the isotopic ratios observed in the nova 
presolar grain candidates. Improving the reaction rates can 
also constrain nova models and simulations, and amend our 
understanding of nova nucleosynthesis [4]. 

According to hydrodynamic classical nova simulations [7], 
the dominant nova nucleosynthetic path is sensitive to the 
chemical composition of the white dwarf, the extent to which 
convective mixing occurs between the material of the white 
dwarf's core and that of the envelope, and the thermal history 
of the envelope. Such questions can be partially answered 
via analysis of the Si isotopic abundance ratios (^^Si/^**Si and 
^"Si/^^Si) in SiC presolar grains of potential nova origin [7], 
and thus such ratios are of specific significance to this work's 
motivation. 

To explore and improve the silicon isotopic abundances in 
nova ejecta predicted from nova simulations, the thermonu- 
clear reactions that most strongly affect the synthesis of sili- 
con in novae must be determined and their rates understood. 
One such reaction is ^^P(/:',y)^''S. Over the temperature range 
characteristic of explosive nucleosynthesis in novae (0. 1 - 0.4 
GK), the rate of the ^^P(/?,Y)^''S reaction competes with that 
of '^P(P+) decay. If in this temperature range the ^'^P(p,j)^^S 
reaction rate is faster than the ^^P(P+) decay rate, and if the 
^''P(P+) decay rate competes favorably with the rate of proton 
capture on ^"P [8], the net effect is an increase in the produc- 
tion of ^''Si via the 2''p(;,,y)30s(p+)30p(p+)30sj reaction se- 
quence, as well as a simultaneous decrease in the abundance 
of ^^Si, which is the product of the p+-decay of ^^P. There- 
fore, an excess in ^"Si together with the depletion in "^Si ob- 
served in some SiC presolar grains could indicate imprints of 
a nova origin. In a study on the sensitivity of nova nucleosyn- 
thesis to uncertainties in thermonuclear reaction rates [1], a 
change in the ^'^P{p,"f)^'^S rate by 10'*, which was consistent 
with the rate limits from Ref. [9], resulted in chang es in 2930si 
abundances by a factor of 3. 

In the temperature range characteristic of explosive hydro- 
gen burning (0.1 < T < 1.3 GK), the Gamow window of the 
2^P(/9 , y)^"S reaction spans £„„ « 700 - 1 770 ke V, where there 
is a low level density. Thus, the rate depends on the properties 
of isolated and narrow ^^P + p resonances corresponding to 
(fi/2 = 1175.9(17) ms [10]) proton unbound states with 
4.5 < £v < 6 MeV. 

The 29p(p,Y)30s rate was evaluated by Wiescher and 
Gorres [11], and more recently by Iliadis et al. [9, 12] and 



Bardayan et al. [13]. The rate calculated by Iliadis et al. [9] 
was found to be dominated by the 3 and 2 J proton unbound 
states in ^^s. The excitation energies corresponding to these 
two unobserved resonances were predicted [9] using the Iso- 
baric Multiplet Mass Equation (IMME) to be 4733(40) keV 
and 4888(40) keV for the states with 7^ = 3+ and 7" = 2+, 
respectively. Such large uncertainties in the resonance ener- 
gies, , resulted in an uncertainty in the rate which spanned 
~4 orders of magnitude [9]. Prior to this prediction, sev- 
eral experiments had been performed to study the structure 
of [14-18]. However, the two astrophysically important 
states predicted by Iliadis et al. [9] were not observed in any 
of the previous experiments. 

A direct measurement of the ^^P(/:',y)30s reaction is cur- 
rently not feasible because no ^^P radioactive ion beam with 
the required beam intensity (> 10** pps) is available. Thus, 
following the prediction by Iliadis et al. [9], attempts were 
made to find these two states via indirect methods [13, 19- 
22]. Bardayan et al. [13] remeasured the excitation energies 
and spin-parity assignments of the states of ^"S up to 7.1 MeV 
by means of the ^^S{p,t)^^^S two-nucleon transfer reaction. As 
a result, a state at 4704(5) keV was discovered and was pro- 
posed to be the predicted 3 state. However, no trace of the 
other important level was found. 

Shortly thereafter, we performed two separate experiments, 
each with two phases, to determine the excitation energies and 
spin-parity assignments of several states of ^''s, which were 
populated via the '^'^S{p,tf'^S and ^^Si(^He, nfr'^S two nu- 
cleon transfer reactions. 

In Ref. [23], the resonance energies corresponding to six 
proton unbound states with < 5.5 MeV in were pre- 
sented, including both astrophysically important states pre- 
dicted by Iliadis et al. [9] one of which was observed for the 
first time. Since then, we have performed a new ^^S(p,t)^'^S 
measurement with a different target (phase II), and have im- 
proved upon the analysis of the existing data. Phase I of our 
^**Si(^He, nyj^'^S experiment was performed with the sole pur- 
pose of determining via Y-ray coincidence measurements the 
energies of the two important resonances predicted by Iliadis 
et al. [9], and phase II was carried out to measure the Y-ray an- 
gular distributions and Y-Y angular correlations from oriented 
nuclei to infer the spins of the observed ^^S states. The results 
of phase I of our ^^Si(^He, nyj^'^S experiment are also pub- 
Ushed [24]. 

The present work discusses in detail the experimental 
setups and data analyses for the second phases of our 
^^S(/7,f)^°S and ^^Si(^He, nyj^'^S experiments, and presents 
our unpublished data for the first phase of our ^^S{p,t)^^S 
experiment. This work thus presents our combined final re- 
sults on the energies and spin-parity assignments of the ob- 
served states, the most updated ^'^P{p,"f)^'^S reaction rate 
calculated via a newly developed Monte Carlo method, as 
well as the impact of this rate on the abundance ratios of 
silicon isotopes synthesized in novae. Therefore, the results 
in the present paper supersede those of our previous publica- 
tions [23, 24]. 
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II. EXPERIMENTS 
A. The ^^S(p, f)^°S experiment 

1. Experimental setup and data analysis: phase II 

The experiment was performed at the Wright Nuclear 
Structure Laboratory (WNSL) at Yale University. A proton 
beam was accelerated, using the ESTU tandem Van de Graaff 
accelerator, to 34.5 MeV (A£/£ - 6 x 10"'*) [23, 25]. 

The beam impinged on a 55.9 ± 5.6 /Jg/cm^ isotopically 
pure (99.9% enriched) '^C foil implanted with 10.4 ± 0.4 
/jg/cm^ of ^^S. This target was fabricated specifically to re- 
duce the relatively flat background produced by the ""'Cd, 
where nat refers to natural, component of the CdS target 
used in phase I of our ^^S(/:',f)^''S experiment [23]. The pro- 
duction procedure for the implanted target is described else- 
where [25, 26]. The thicknesses of the ^^S and '^C layers 
in the implanted target were obtained through a Rutherford 
backscattering measurement [25, 26]. 

In addition to the aforementioned target, a free-standing 
311-/jg/cm^ natural Si foil was used for calibration purposes. 
Also, a stand alone 40-/jg/cm^-thick 99.9% isotopically en- 
riched '^C foil was used to measure the background from (p, f ) 
reactions on the carbon substrate in the implanted target. The 
method of measuring the thicknesses of these targets is de- 
scribed in Ref. [25]. 

The reaction ejectiles were dispersed according to their mo- 
menta with an Enge split-pole magnetic spectrograph, with 
vertical and horizontal aperture settings of A(|) = ±40 mrad, 
and A9 = ±30 mrad, respectively. The study was carried out 
at multiple angles with magnetic field strengths of 10 kG for 
9 = 22°; 9.5 kG for 9 = 27.5°; and 9.2 kG for 9 = 45°, where 
9 is the scattering angle in the laboratory system. 

The tritons were focused at the spectrograph's focal plane, 
where they were detected with an isobutane-filled position 
sensitive ionization drift chamber [25], together with a plastic 
scintillator. The ionization chamber measured the positions 
along the focal plane and energy losses {AE) of the tritons. 
Those that passed through this detector deposited their resid- 
ual energy (iSm) in the plastic scintillator. 

A/i, Eyes and position (proportional to momentum) were 
measured to identify tritons and determine their momenta. 
The tritons were selected according to AE and Eres, which 
were plotted vs. focal plane position gates. The spectra of the 
tritons' momenta were then plotted for each spectrograph an- 
gle (see Fig. 1). Triton peaks corresponding to ^"S states in 
these spectra were clearly identified through kinematic analy- 
sis. 

The only contaminant peak observed was the ground state 
of ^'^C (see Fig. 1). The ^'S implanted target produced a back- 
ground that was decreased by about a factor of 2 compared 
with the relatively flat background produced by the ""'Cd com- 
ponent of the CdS target used in phase I of the ^^S{p ^t)^^S 
experiment. 

The triton peaks observed in the presented spectra were fit- 
ted using a least-squares multi-Gaussian fit function to deter- 
mine the peak centroids, widths and areas. The energy cali- 




6 7 8 9 10 

Triton Lab Energy (MeV) 



Figure 1: Triton spectra from the S{p,t) S reaction obtained with 
the implanted target. Peaks corresponding to ■'''S states are labeled 
with energies in keV. The filled histograms are background spectra 
measured with an isotopically enriched '^C target, normalized to the 
"^^S(p,0^'^S data. The main contaminant is the ground state (g.s.) of 
"'C. For 27.5°, an aluminum plate along the focal plane blocked the 
region corresponding to tritons with energies higher than 9.5 MeV, 
where elastically scattered protons reached the focal plane. At 45° 
the gates cut the region to the right of the peak corresponding to the 
4688-keV state. 

bration was determined from a combination of known levels 
of ^^Si (measured with the ^^Si(/:',f) reaction using the Si- 
target) and of ^"S, whose adopted energies are weighted av- 
erages of previous work on ^^Si levels [14, 27-33] and on 
^"S [14-17]. Since the earlier publication [23], the previous 
calibration fits were improved through reanalysis of the pre- 
vious data (corresponding to phase I) by accounting for the 
angle of the target with respect to the beam (details are pro- 
vided in Ref. [25]). Fig. 2 presents the Ex > 5.5 MeV excited 
states in '^"S observed in phase I of the ^^S{p,t)^'^S experiment 
that were not published in Ref. [23]. 

The final excitation energy uncertainties for the data of both 
phases of this experiment arise from: (1) statistical uncer- 
tainties (< 2 keV), (2) uncertainties in the thicknesses of the 
implanted target (1 keV) and the CdS target (2 keV) taking 
into account the uncertainty in the thickness of the ""'Si target 
used for calibration, and (3) uncertainty in the Q-values of the 
^^Si(p,t)^^Si and ^^S{p,t)^'^S reactions (0.3 keV [34] and 0.4 
keV [10], respectively). Therefore, the ^"S excitation energy 
uncertainties, when added in quadrature, were 3 keV and 2 
keV for the CdS and implanted targets, respectively. Lastly, 
to obtain the final ^*^S excitation energies, a weighted average 
was calculated for each state over all the angles, and thus over 
both targets. 

With respect to the previous publication [23], all the mea- 
sured excitation energies from the present work have smaller 
uncertainties by at least 40% as a result of a reduction in the 
uncertainty of the Q-va\ue of the ^^S{p,t)^'^S reaction due to a 
recent improved measurement [10] on the ^"S mass, and our 
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Figure 2: Triton spectra measured (in phase I) from the S(p,t) S 
reaction at selected angles obtained with the CdS target (for details, 
see Ref . [23]). Peaks corresponding to ■'"S states are labeled with en- 
ergies in keV. The states with £■( > 5.5 MeV were not published 
in Ref. [23]. The filled histograms are background spectra mea- 
sured with a ""'Cd target on a carbon backing, normalized to the 
^^S(p,0""'S data. Apeak from the '■'C(p,0"C reaction is also iden- 
tified and labeled by its parent nucleus. For 10° and 20°, an alu- 
minum plate along the focal plane blocked the region corresponding 
to tritons whose energies are higher than 1 1 MeV, where elastically 
scattered protons reached the focal plane. 



improved calibration fits for the previous data obtained by the 
CdS target. 

The energy resolution was approximately 28 keV and 22 
keV (FWHM) for the spectra obtained with the CdS and im- 
planted targets, respectively. Therefore, our achieved energy 
resolution is a factor of 3 - 5 smaller than those of previous 
^^S(p,tf^S measurements [13, 14]. 

2. Results: both phases combined 

Over both phases of the ^^S{p,t)^^S experiment, 12 proton 
unbound states of ^"S with < 6.8 MeV were observed, and 



Table I: Weighted average excitation energies of S from both 
phases of our ■'^S(p,f)""'S experiment. States used for energy cali- 
bration are marked by an asterisk. 



Ex (keV) 


J"" 


Ex (keV) 




2208(3) 




5393(2) 


3+ 


3402.6* 




5849(2) 


(1^,2+, 4+) 


3681(3) 


(l+,0+) 


5947(2) 




4688(2) 


3+ 


6055(3) 


(1-) 


4812(2) 


2+ 


6345(3) 


(0+) 


5136* 


(4+) 


6536(3) 


(2, 3) 


5225(2) 


(0+) 


6768(3) 


2(-) 


5315(2) 


(3-, 2+) 







their weighted average energies (over all angles) are listed in 
Table I. 

Most of the measured energies in the present work are in 
agreement within 1 - 2o with those measured in the previ- 
ous ^2s(;j,f)30s measurements [13, 14]. The 5947-keV state 
observed in the present work (see Fig. 2) is in good agree- 
ment with the 5945-keV tentative level observed in Ref. [18]. 
For the state with an expected excitation energy of « 4.7 
MeV [9], our measured energy of 4688(2) keV does not agree 
with the 4704(5) keV measured in Ref. [13]. Most of the lev- 
els observed in our ^^S{p,t)^'^S experiments whose > 5 
MeV have been measured previously but have spin-parity as- 
signments that are either unknown or tentative. 

To obtain the spin-parity assignments of ^"^S states observed 
in both phases of the ^^S{p,t)^'^S experiment, the equiva- 
lent thickness of the sulfur content of the CdS target was re- 
quired. This thickness was determined to be 53 ± 5 /Jg/cm^ 
through the reanalysis of the data of a previous scattering ex- 
periment [35], where an 8-MeV '*He+ beam along with the 
Enge spectrograph at WNSL and a silicon surface barrier de- 
tector were used to determine the composition and thickness 
of the CdS target. The theoretical angular distributions of 
the cross sections were then computed, following the work 
of Ref. [13], using (i) the one-step finite range transfer for- 
malism via DWUCK5 [36] for the natural-parity states, and (ii) 
the coupled reaction channels formalism via FRESCO [37] (un- 
der the assumption of finite-range interaction potential) for the 
unnatural-parity states. DWUCK5 is unable to handle unnatural- 
parity states. FRESCO, on the other hand, can be used for all 
states. The former code was used for the cases of natural- 
parity states due to its simplicity. Its output for a few natural- 
parity states was examined against that of FRESCO, and the 
results were identical in terms of the shapes of the theoretical 
DWBA curves. 

To determine the orbital angular momenta that were trans- 
ferred in the ^^S(/7,f)^*'S reaction to populate the observed 
^"^S excited states, the same shell-model assumptions made in 
Ref. [13] are made here. The angular distributions for the se- 
quence ^^S{p,d)^^Sf,,sXd ,t)^'^S were obtained for each of the 
unnatural -parity final states in ^"S. 

The distorted waves in the entrance and exit channels were 
calculated for optical interaction potentials, the parameters of 
which were taken from Ref. [13] (and references therein), and 
are given in Table II. 
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Table II: Optical model parameters used in conjunction with DWUCK5 and FRESCO for the analysis of the angular distributions. 



Reaction Vq Wq Wd vq a r'^ a' r(( a" "k PNLOC 

Channel (MeV) (MeV) (MeV) (MeV) (fm) (fm) (fm) (fm) (fm) (fm) (fm) (fm) 





37.1 





6.875 


7.5 


1.18 


0.66 


1.18 


0.66 


1.18 


0.7 


1.25 




f + 30s 


144 


30 








1.24 


0.68 


1.45 


0.84 








1.25 




rf + 3is 


90 





25 




1.30 


0.62 


1.18 


0.58 






1.25 




n + 3is'' 












1.20 


0.65 










1.30 


25 


2n + 












1.25 


0.65 












25 



"ThLs parameter, a.s well as and a" , is taken from Ref. [14]. 

'The input parameters corresponding to this channel are taken from Ref. [38]. 



Furthermore, the widely used Reid soft core potential [39] 
was used to derive the deuteron and triton wave functions, as 
well as the p-n and d-n interactions for FRESCO calculations. 

Triton angular distribution plots are shown in Fig. 3. An- 
gular distributions of the states with £, < 5.136 MeV are dis- 
cussed in § IIB2(2.4), since those states were also observed 
in our y-ray measurements. In the following, we will only dis- 
cuss the angular distributions for levels with > 5.225 
MeV. 

• The 5225-keV level: This state is a prominent peak that 
was observed at every angle measured in the ^^Si^p^t^'^S ex- 
periments. There is no conclusive information regarding the 

assignment of this state in the literature. Our only guide 
comes from a shell-model calculation [11], which suggested 
that there should be a 0+ level around 5.2 MeV. Although rea- 
sonable fits are obtained with /"^ = 0+ and 2+, our fit using 

= 0+ is more consistent with the data. So we suggest 
that this state is the mirror to the state in ^^Si at 5372.2 
keV [40]. 

• The 5315-keV level: This state is also a prominent peak 
observed at all angles. It is known to be a 3^ state [17]. Our 
angular distribution is better fitted by an Z = 3 angular mo- 
mentum transfer, but I = 2 would also be reasonably consis- 
tent (see Fig. 3). We suggest that this state is most likely the 
mirror to the 5487.5-keV state in ^^Si with 7" 3" [40]. 

• The 5 39 3 -keV level: This state was observed at all angles 
measured in the ^^S{p^t)^^S experiments. Its spin was ten- 
tatively assigned to be 7 = 1 or 2 in previous work [17]. In 
Ref. [13], tentative = 3^ and 2+ assignments were given 
to this state. Our triton angular distribution is more consis- 
tent with = 3+ assignment, and thus we assign this state to 
be 3+, making it the mirror to the 32^ state in ^''Si at 5231.38 
keV [40]. 

• The 5849-keV level: This state was tentatively assigned 
to be a 1" state in Ref. [13]. However, / = 2, 3 and 4 transfers 
could not be excluded. In our data, this level was observed at 
10°, 20°, 22° and 45°. We can rule out 7" = 4" and 2" as- 
signments but 1^, 2+ and 4+ are all in reasonable agreement 
with our data (see Fig. 3). 

• The 5947-keV level: This level was too weakly populated 
to extract a significant angular distribution. 

• The > 6 MeV states: With the exception of the 6055- 
keV and 6768-keV states, which are observed at four angles 



(see Fig. 3), all other states of S observed in the present 
experiment whose excitation energies are above 6 MeV are 
only observed at most at three angles, 10°, 20° and 22° (see 
Fig. 2). Nonetheless, we propose a tentative assignment of 
1^ to the 6055-keV state, which is consistent with the results 
of Ref. [17], but the energy of this state from our data differs 
by 62 keV. We also make a tentative assignment of 0+ to the 
6345-keV state consistent with a definite = 0+ assignment 
made in Ref. [17]. Furthermore, we tentatively assign J — 2 
or 3 to the 6536-keV state, which is consistent with what was 
suggested in Ref. [ 17] . Lastly, for the 6768-keV state, we con- 
firm 7 = 2 suggested in Ref. [13], and likely rule out / = 3 and 
4 transfers. Though, our angular distribution data are best fit- 
ted with a negative parity assignment. Therefore, we propose 
a 7" = 2^^' assignment to this state. 



B. The 2^Si(3He,«Y)30s experiment 

1. Experimental setup and data analysis: phase II 

An in-beam y-ray spectroscopy experiment using the 
^**Si(^He, «y)-"'S reaction was carried out to assign spins to 
the populated ^^S levels based on measurements of y-ray an- 
gular distributions and y-y angular correlations from oriented 
nuclei. This experiment was performed at the University of 
Tsukuba Tandem Accelerator Complex (UTTAC) in Japan. A 
^He^+ beam was accelerated to 9 MeV via the 12UD Pelletron 
tandem accelerator at UTTAC. The details of this beam are de- 
scribed in Ref. [24]. The beam impinged on a self-standing 25 
^m-thick foil of ""'Si, of which the ^^^Si abundance is 92.23%. 

High-purity germanium detectors with 50% and 70% rel- 
ative efficiency were placed at 90° and 135° with respect to 
the beam axis, respectively. We hereafter refer to these de- 
tectors as 1 and 2, respectively. These detectors were located 
on opposite sides with respect to the beam line. The energy 
resolution of detectors 1 and 2 was determined to be 4.4 keV 
and 3.2 keV (FWHM) at Ey = 1333 keV, respectively, y-y co- 
incidence data were accumulated during a total of 4 days, and 
was corrected hourly for detector gain shifts. A sample y-ray 
spectra can be seen in Fig. 4. 

To extract the centroid and area of each peak, the peaks 
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Figure 3: Triton angular distributions populating states of S compared with the DWBA curves obtained by DWUCK5 and FRESCO for natural- 
and unnatural-parity levels, respectively. The filled circles with error bars are the measured differential cross sections in the center-of-mass 
system, and the solid, dashed or dotted curves are the theoretical angular distributions. If not shown, the error bar is smaller than the point size. 
The excitation energies are given on the top middle of each plot. The shell-model orbital from which the di-neutron is stripped is given for the 
0+ and the negative natural-parity levels. 
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(a) y-y coincidence spectrum measured at 90° 



Figure 4: Singles y-ray spectrum measured during pliase II of the 
experiment at 90° using detector 1. Selected strong transitions are 
labeled by their parent nuclei and with energies (in keV) that are 
weighted averages between both phases of the experiment. The 
2210.6-keV (2;^ O'l) and 1194-keV (2 J 2+) peaks originate 
from levels in -"'S. 



were fitted using a single-Gaussian function whenever they 
were reasonably isolated from each other, and with a multi- 
Gaussian function for the partially resolved or unresolved 
doublets. Those peaks that were affected by Doppler shift 
at higher angles were fitted using Gaussian-plus-exponential 
functions to account for the exponential tail. Background sub- 
traction was performed by assuming a linear function under 
each peak. 

The Ge-detectors' initial energy calibration and energy- 
dependent efficiencies were determined with a standard ^^^Eu 
calibration source. The initial energy calibration fit was im- 
proved via internal calibration by using strong ^"P y-rays emit- 
ted from the ^**Si(-'He, py)-'"? reaction, whose cross section is 
higher than that of the ^^Si(^He, «y)^*^S reaction at this beam 
energy [41, 42]. The resulting uncertainties in the detection 
efficiencies were estimated to be 5%. 

The coincidence analysis was performed via construction 
of a y-y coincidence matrix. Fig. 5 presents the coincidence 
spectra. The y-ray angular distribution and y-y angular cor- 
relation measurements and their results will be discussed in 
§ II B 2 (2.2) and § II B 2 (2.3), respectively. 



2. Results: both phases combined 

(2.1). Decay scheme of ^^S 

In the singles y-ray spectra of both Ge-detectors during each 
phase of the experiment, two y-rays were clearly observed at 
2210.6(3) keV and 1194.0(1) keV, which correspond to the 



2+ 



0^ 



and 2J — s- 2^ transitions in -^''S, respectively (see 



Fig. 4). 

A few y-rays with energies in the range of 3 MeV were ex- 
pected to be observed in the singles spectra according to the 
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(b) y-y coincidence spectrum measured at 135° 

Figure 5: The y-y coincidence spectrum measured during phase II 
at (top) 90° and (bottom) 135° obtained from gating on the 2210.6- 
keV (2j^ — Oj^) transition of -"'S. Peaks corresponding to the tran- 
sitions from known ■''^S states are labeled with energies (in keV). At 
135°, the labeled energies are corrected for Doppler shift except that 
of the 846.1-keV y-ray (see text). The 2477.1-keV and 2599.5-keV 
peaks are from the decays of proton-unbound states at 4688.0 keV 
and 4810.4 keV, respectively. 



measured branching ratios [15, 16] of the y-rays from decays 
of the bound states and the lowest-lying resonances of ^''S. 
However, these y-rays did not appear as separate observable 
peaks in the singles y-ray spectra obtained during either phase 
of the experiment. This was most likely because they were ob- 
scured by the Compton scattered y-rays from ^'"P transitions. 

After placing software gates on the 2210.6- and 1194.0- 
keV peaks, y-decay cascades from higher-lying states were 
observed in the y-y coincidence spectra (see Fig. 5). In partic- 
ular, we observed transitions with energies of 2477.3(3) keV, 
2599.0(4) keV (see Fig. 2 in Ref. [24]) and 2921.4(4) keV 
from ^^S proton-unbound states at 4688.1(4) keV, 4809.8(5) 
keV and 5132.3(5) keV, respectively. 

Recoil energies were taken into account when construct- 
ing the final excitation energies of ^"S from its y-ray decay 
scheme. The results are given in Table UI. The final uncer- 
tainties in the energies are due to the statistical uncertainties 
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Table III: Weighted average energies (between both phases of the experiment) and relative intensities of the observed transitions in S. The 
latter are calculated with respect to the strongest y-ray measured at the same angle. The uncertainties in the recoil energies (E,.^^„ji) were 
negligible, and thus are not presented. The energies of initial and final states {E, and Ef, respectively) are corrected for the corresponding 
recoil energies. The results obtained in the y-ray measurement of Ref. [16] are also shown for comparison. 



Present Work Ref. [16] 
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^recoil 
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ly 
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Ei 






(keV) 


(keV) 


(keV) 


(keV) 


(%) 


(%) 


(keV) 


(keV) 


(keV) 
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846.0(4) 


0.01 


unplaced 
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3.9(6) 


2.8(5) 










1194.0(1) 


0.03 


3404.7(3) 
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33.5(5) 


43.3(10) 
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2210.7(5) 
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0.03 
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0.04 


3667.2(4) 


2210.7(3) 


11(3) 


13.9(9) 


1456.8(9) 


3667.5(10) 


2210.7(5) 


2210.6(3) 


1466.2(3) 
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0.10 


4688.1(4) 


2210.7(3) 
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0.10 
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"The transition on which the coincidence gate is placed. 
'This transition is not observed at 135°. 

''This transition is too weak at 135° to obtain a reasonable yield. 
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Figure 6: (left) A portion of Si decay scheme [40] in comparison with (right) that of its mirror nucleus, S, based on the results of the 
present work. The y-ray branches are not to scale; however, the thicker the arrow, the stronger the branch. All the observed y-rays are shown 
with their energies (in keV) corrected for the recoil energies of the corresponding -"'S excited states, which are also shown (in keV). The y-ray 
transitions with energies above 3 MeV in ■'"S could not be resolved in our experiments due to the presence of strong ■'''P transitions in that 
region. 



in the corresponding centroids only, because all the calibra- 
tion energies have negligible uncertainties. 

From the recoil energies and the y-ray energies, the exci- 
tation energies of the first few states were reconstructed to 
obtain the level scheme of ^"S (see Table III and Fig. 6). The 
measured energies of most of the levels are in agreement with 
the results of the ^^S(p,t)^^S measurements discussed earlier, 
as well as those of previous measurements on the y-rays of 
^"S [15, 16]. In particular, the measured energies of the two 
astrophysically important excited states at 4688.1(4) keV and 
4809.8(5) keV from our y-ray measurements are in excellent 



agreement with the 4688(2) keV and 4812(2) keV energies 
from the ^^S(p,t)^'^S experiments presented in Table I. 

However, there are discrepancies in the energies of two ^"S 
levels: the energies of the 3404.7- and 5132.3-keV states de- 
duced from our y-ray energies are ~2 keV higher and ~4 keV 
lower, respectively, than those measured in the y-ray study of 
Ref. [16]. The reason for the discrepancy in the energy of the 
5-MeV state is unclear; however, we suggest that the inconsis- 
tency between the measured energies of the 3-Me'V state orig- 
inates from the presence of a double escape peak at 1188.6 
keV (see Fig. 7) just beside the peak at 1 194 keV, correspond- 
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Figure 7: (Color online) The 1194-keV peak together with the dou- 
ble escape peak of the 2210.6-keV line of ^''S. The other two peaks 
are identified by their parent nucleus and their energy (in keV). The 
black and red spectra are measured by the 50% and 70% relative effi- 
ciency detectors, respectively. For a short time during Phase II of the 
experiment, these detectors were placed at ±90° with respect to the 
beam axis. The I188.6-keV transition is the double escape peak of 
the 2210.6-keV y-ray of -"'S, and its yield has decreased significantly 
when measured by the larger detector at -90°. The peak correspond- 
ing to the 1194-keV y-ray is one of the two that stands out in the 
spectrum measured by the larger detector at -90°. 



ing to the 2^ ^ 2j transition in ^'^S, observed in our sin- 
gles spectra. 

The energy of the 1194-keV y-ray results in the level en- 
ergy of the 2 J state of ^"S to be ~2 keV higher than that 
measured by Kuhlmann et al. [16]. The latter measurement 
was performed in the early 1970's when the Ge-detectors were 
smaller Thus, it may be possible that the 1 188.6-keV double 
escape peak was also present in their spectra but because of 
the lower detector efficiency the two peaks were assumed to 
be one. 

We expected to observe the y-rays emitted from de- 
excitations of the 3407.7- and 3676.9-keV states directly to 
the ground state in the singles spectra. In addition, if the 
4809.8-keV state is the 2\ state in ^^S, then according to the 
decay scheme of its mirror level we expect that the 2^ — > Oj^ 
transition in ^"S is a strong branch (with respect to the strength 
of the other decay branches of the same level). Therefore, we 
also expected to observe the 4809.8-keV y-rays of ^"S in the 
singles spectra. However, the detection efficiency for detect- 
ing such high energy y-rays is relatively low, and the high en- 
ergy regions of the spectra obtained in the ^^Si('^He, «y)^''S 
experiments are obscured mostly by wide peaks originating 
from transitions in ^*^P. Therefore, the 3407. 7-, 3676.9-, and 
4809.8-keV transitions are not resolved. Hence, the fact that 
the 2^ — > 0[^ transition in '^"S is not observed in these exper- 
iments does not imply that this transition is weak. Based on 
Ref. [43] where the intensities of the y-rays of the mirror nu- 
cleus -"^Si were measured at 90°, we estimate that the branch- 
ing ratio of the 2^ — > Oj^ transition in ■"'Si is 36 ± 3 % which 
should be similar to that of the transition from the 4809.8-keV 



state to the ground state in ^"S. 

We have observed a weak line at 846 keV in the singles 
y-ray spectrum measured at 90° (see Fig. 4), which also ap- 
pears in the coincidence spectra at 90° and 135° as a more 
noticeable peak (see Fig. 5). The energy of this peak does not 
seem to be Doppler shifted at 135°, which suggests that this 
y-ray may originate from a state whose half-life is more than 
2 ps [44]. This y-ray transition is also in coincidence with the 
1194-keV transition in ^"S. A weighted average between in- 
dependent measured energies at 90° and 135° for this y-ray 
results in Ey = 846.0(4) keV (see Table III). The energy of 
this transition does not add up to any of the known levels of 
^"S; however, the fact that it is a fairly prominent peak and is 
in coincidence with two transitions of ^''S suggests that this 
y-ray may also belong to the decay scheme of this nucleus. 
The higher lying resonances {E^ > 6 MeV) of ^"S may emit 
y-rays in this energy range, e.g., the decay transition from the 
state with E;, 7123(10) keV to that with E,r = 6280.1(12) 
keV [40]. In particular, if one of these resonances has a high 
spin, its proton decay might be suppressed by the centrifugal 
barrier, and thus it can decay via y-ray emission. The 846- 
keV y-ray transition has so far remained unplaced in the level 
scheme obtained from the present experiment. 

The relative intensities from full-energy-peaks of all the ob- 
served transitions were calculated at 90° and 135°. For the 
coincidence spectra, first the yield of the 1 194-keV transition 
observed in the singles spectrum was normalized to that of 
the 2210.6-keV y-ray transition also obtained from the same 
spectrum. Then the relative intensity of the 1194-keV y-ray 
transition was used to convert the yields of all the other y-ray 
transitions in the coincidence spectra into relative intensities. 
These results are tabulated in Table III. 

(2.2). Singles measurements: y-my angular distributions 

For a transition 7, Jf, where J represents the spin of the 
state, the theoretical y-ray angular distribution function is de- 
fined as [45]: 



W{Q)rheo = E AkPkicosQ), 

k = 
k = even 



(1) 



where the coefficients 7\ (cos9) are the Legendre Polynomi- 
als, and the coefficients are defined as [46] (p. 55): 

+ 5^Fk{jfXrj,)], (2) 

where y, and jf are the spins of the initial and final states in- 
volved in the transition, respectively; A, and A,' are transition 
multipolarities; a^; are the alignment factors (see Eq. (4)); 
and Fk coefficients are tabulated [47] for different y, — s- jf 
transitions; and 5 is the mixing ratio of a y-ray transition de- 
fined as [48] (p. 54): 
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Figure 8: Experimental y-ray angular distributions of the 2210.6-keV (left) and the 1194-keV (right) transitions. Both y-rays are observed in 
the singles spectra obtained during phase 11 of the experiment. They correspond to the ^ ""^^ ^ trQr,c;t;^r,c in 30i; 
The solid lines are best fits to Legendre polynomials. 



■ Oj and 22 — >■ 2j transitions in S, respectively. 



where co'L and coL + 1 are the magnetic and electric transi- 
tions of multipolarity L, and L + I, respectively; and Fj is the 
partial y-ray width corresponding to a state with spin j. 
The aUgnment factors are defined as [46] (p. 55): 



(4) 



where k is even and k > 6 are ignored due to a rapid de- 
crease of transition probabilities of higher order multipoles. 
An individual aligned state with spin j can be represented as 
a Gaussian probability distribution P{m) of 2^ + 1 magnetic 
substates nij, where mj = —j, • • • , j, with the FWHM of o 
along the beam axis. P{m) is the population parameter and is 
defined as [46] (p. 56); 



exp 



P{m) 



2o2 




(5) 



where o can be defined experimentally. 

An incomplete alignment of a state relative to the beam 
axis results in an attenuation of the population parameter The 
alignment factors, a^ , describe the degree of the attenuation of 
the population parameter. While ao is considered to be unity, 
0,2 and a4 coefficients are determined experimentally. 

For the y-ray angular distribution measurement during 
phase II of the experiment, the total charge deposited by the 
beam could not be determined due to a faulty beam current 
integrator. Therefore, to take into account the fluctuations in 
the beam intensity and possible target degradations or changes 
in the target profile that could affect the areas under the peaks 
of interest, detector 1 was used as a monitor detector. It was 



kept fixed at 90° with respect to the beam axis 10 cm away 
from the target. Detector 2, on the other hand, was positioned 
7 cm away from the target and on the opposite side of detector 
1. Detector 2 was moved between 90° to 120° in intervals of 
10° and was lastly positioned at 135° with respect to the beam 
axis. It could not be place at angles higher than 135° due to 
the presence of the beam line. 

The singles y-ray spectra were then obtained for 1 hour 
from both detectors at five different angular pairs and were 
calibrated as explained before. 

For every (9i, 82) angular pair corresponding to detectors 
1 and 2, the intensities of the 2210.6-keV and 1 194-keV y-ray 
transitions from ^"S were normalized to the intense 1+ — > 1 + 
transition at 708.7 keV in ^°P. 

The normalized relative yields for each peak of interest 
were plotted against cos2(9), where 9 is the detection angle, 
and these data were fitted (see Fig. 8) using the function: 



W{Q)e;,p = Aq + A2P2{cosQ) +A4P4(cos9), 



(6) 



where the coefficients A, are extracted from the fit, and 
P2(cos9) and P4(cos9) are Legendre polynomials. W{Q)exp 
represents the experimental y-ray angular distribution func- 
tion, which can be used to normalize W{Q)theo- From the lat- 
ter, one can infer the alignment probability of an excited state 
involved in a y-ray transition. 

The angular distributions of the two observed ^"S peaks in 
the singles spectra are discussed below. 

• The 22 10.6-keV transition of^'^S: 

The 2210.6-keV y-ray corresponds to the 2'^ — >■ 0j~ tran- 
sition in ^"S. This transition is a pure E2 (8 = 0) and is a 
so-called stretched quadrupole transition [49]. 

The experimental intensities of the 2'^ 0^ transition (see 
Fig. 8) was used to normalize the W{Q)theo of this transition 



11 



Table IV: Results of the y-ray angular distribution studies for -^"S transitions observed in the present work. Energies are in keV. 
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0.38(25) 
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E2 

Ml, £2 




0.16 



"This value is normalized such that Eq. (6) becomes W„p(9) = 1 + {A2 / Ag) Pi^cosQ) + (A4//lo)P4(cos9), which resembles Eq. (1), where AoPo(cos9) = 1. 

*The attenuation factors are [50] (a2,a4) = (0.41482,0.048393) for a/j = 0.6, and (a2,a4) = (0.53784,0.095181) for c/j = 0.5. 
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Figure 9: Experimental y-ray angular distributions shown in circles in comparison with the theoretical angular distributions normalized to the 
data shown with solid lines. The former were obtained from fitting Eq. (6) to relative intensities shown in Fig. 8, and the latter were calculated 
using Eqs. (1) and (2). The comparison is made for the (left) 2210.6-keV and (right) 1194-keV y-rays. The normalization of the theoretical 
angular distributions is best performed with ct/ j = 0.6 and 5 = for the 2210.6-keV y-ray, and with 0/ j = 0.5 and 5 = 0. 16 for the 1 194-keV 
y-ray. The agreement between the theoretical and experimental curves over most of the angles is good for the 2210.6-keV y-ray and less 
satisfactory for the 1 194-keV y-ray (see text). For those angles lower than 40° and higher than 150°, the VV(6)™™' diverges significantly from 
W(Q)exp due to the lack of data points for normalization at those angles. 



obtained via Eqs. (1) and (2) using 8 = 0. 

To normalize W{Q)theo to W{Q)exp, the coefficients B2F2 
and B4F4 for the 2+ 0+ transition were taken to be 0.7143 
and -1.7143, respectively, from Ref. [46] (p. 82). Therefore, 
the only parameters that were free to vary were the alignment 
factors a2 and a4. These coefficients are given in the litera- 
ture [50] for 0.1 <a/j < 2. 

Thus, for each (a2, a4) pair corresponding to a specific a/j 
value, the theoretical angular distribution was calculated at the 
same angles at which a relative yield was measured in phase II 
of the ^^Si('^He, «y)^'^S experiment. An average normalization 
factor was thus obtained and was used to normalize W{Q)theo 
to the intensity at each angle. Then \W{Q)fm - hxpWexp, 
where lexp and hhxp are respectively the intensity and its un- 
certainty obtained from the data at the angle 9, was plot- 
ted against cos^ 9. Hence, the specific (a2, OC4) pair, which 
yielded the minimum difference between W{Q)theu and lexp, 
was found. Finding the (a2, pair uniquely determines the 
parameter 0/ j, where o is the FWHM of the population pa- 



rameter The results for the 2^ — !> Oj^ transition in S are 
presented in Table IV and Fig. 9. 

• The 1 194-keV transition of^^S: 

From a comparison of the 1 194-keV y-ray, corresponding 
to the 2^ — !> 2^ transition in ^"^S, with the mirror transition in 
^"Si, it was assumed that this transition is a mixed Ml/£'2. For 
this transition, the mixing ratio 8 is an additional free param- 
eter that is required for normalization of W{Q)theo to W{Q)exp- 

For the 1 194-keV y-ray, the coefficients A2/A0 and A4/A0 
were first extracted from the experimental fit (see Eq. (6)). 
W{Q)rheo was calculated for all (a2, 0C4) pairs corresponding 
to 0.1 < o/j < 2 for a 2+ ^ 2+ transition [50]. The mixing 
ratio was set to a constant free parameter from a prechosen 
set of values. The parameters ^2^2 and B4F4 are constants 
given in the literature [46] (p. 82). With these, W{Q)ti,eo was 
calculated for each value of 8. Then, a statistical test was 
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Figure 10: (Color online) (a) (solid line) vs. arctan6 for the 1 194-keV y-ray de-exciting the 3404.7-keV state of ■''^S. The dashed line shows 
Xm/ii + 1' therefore our la confidence level (see text), (b) and (c) Polynomial fits of the 3*^^ degree are shown with solid red lines passing 
through a selected portion of (circles) vs. arctan6. 



performed with x defined by: 



\ theo 



Mr 



.exp 
^4 



i theo 



(7) 



where the A^^'' parameters are the yields of ^"S y-rays ob- 
served in the singles spectra and normalized to that of a ^'^P 
y-ray peak as discussed earlier; 54^^'^'' are the experimental un- 
certainties in A, normalized to 54o extracted from the fit given 
by Eq. (6); and Af"'" is calculated using Eq. (2). 

The was plotted against arctanS (see panel (a) in Fig. 10) 
and had local minima at arctanS ~ 10 and arctanS ~ 60. 

was again separately plotted for two regions around these 
minima, and each region was fitted (see panels (b) and (c) 
in Fig. 10) with a polynomial of the third degree to obtain 
the functional forms of with respect to 8 for these regions. 
Then a minimization procedure was used to find the best 
possible local solutions, which were 8 = 0.16 and 8 = 1.5. 

Those 8's that are within 1.0 of the best x^ (see the dashed 
line in panel (a) in Fig. 10) are located at approximately ±la. 
Therefore, all 8's within -0.13 < 8 < 3.73 are valid, which 
means our uncertainty in delta is very large. However, our 
choice of 8 from the aforementioned range is determined by 
the consistency with the mixing ratio of the mirror transition 
(8 = 0.18(5) [40]) and the agreement between the theoretical 
and experimental angular distributions for the 1 194-keV y-ray. 
The latter is best for a/ J = 0.4; however, 8 in that case is cal- 
culated to be 0.04, which is not consistent (within 2a) with the 
mixing ratio of the mirror transition. Therefore, the next best 
value is a/ j = 0.5, for which 8 = 0. 16 consistent with that of 
the 2 J — > 2]'" mirror transition in ^"Si. 8 = 0. 16, as mentioned 
before, also represents a local minimum in the vs. arctan8 
plot. We thus adopted 8 = 0. 16 and held it fixed. For the sign 
of 8, we have followed the convention adopted by Krane and 
Steffen [51] as opposed to that of Rose and Brink [52]. 

To confirm that we can reject 8 = 1.5, the single particle E2 
transition strength B(£'2;2+ — > 0+) in Weisskopf units was 
determined as follows: 

^(E2)(inW.u.)= 9.527 xlO^^^ 



where f i /2 is the half-life of the state under consideration; A 
is the mass number; Ey is in keV; a is the internal conversion 
coefficient, which is ignored for our case as this coefficient 
decreases with increasing Ey; and BR is the branching ratio of 
the transition of interest. 

We obtained B{E2;2+ ^ 0+) = 0.41 for the 1 194-keV y- 
ray transition. As a rule of thumb [44], if the B(£'2;2+ ^ 0+) 
of a transition is larger than one, the corresponding state which 
emits the y-ray of interest is most likely a collective state, for 
which the mixing ratio should be large. On the other hand, 
when a transition has B(£'2;2+ — > 0+) < 1, the state which 
initiates the transition is to a good approximation estimated 
as a single particle state with a small mixing ratio. Since our 
estimated B(£'2;2+ — !• 0+) value for the 1 194-keV transition 
falls into the latter category, we concluded that the 3404.7- 
keV state is a single particle state with a small mixing ratio. 
Therefore, we adopted 8 = 0.16. 

Finally, the procedure which was described for the 2210.6- 
keV y-ray transition was repeated for the 1 194-keV y-ray tran- 
sition to determine the FWHM of its population parameter 
The results are given in Table IV and Fig. 9. 

The previously described x^ method was also performed as 
a check for the 2210.6-keV transition, and a sharp minimum 
at 8 = confirmed the stretched E2 profile for this y-ray. 

The large uncertainties in the experimentally determined 
A2/A0 and A4/A0 (see Table IV) are mostly due to the low 
statistics in each peak in the singles y-ray spectra. However, 
they are still consistent with the typical values [44, 53, 54] 
expected for a stretched quadrupole with AJ = 2 (for the tran- 
sition from the 2210.7-keV state to the ground state) and a 
mixed dipole-plus-quadrupole with AJ = (for the transition 
from the 3404.7-keV state to the 2210.7-keV state). 



(2.3). Coincidence measurements: y-y angular correlations 

Measurements of the Directional Correlations of y-rays de- 
exciting Oriented states (DCO ratios) allow to deduce the an- 
gular correlation information from the y-y coincidence data. 
The method of measuring DCO ratios is applied to determine 
the multipolarities of the y-rays involved in a cascade, and thus 
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Table V: The experimental DCO ratios for S y-rays observed in the present experiment, ji, j,„ and jf are the spins of the initial, intermediate 
and final states, respectively. Theoretical DCO ratios are from Refs. [44, 54], and are obtained from known transitions for which a/ j — 0.3. 
See text for further explanations of the theoretical ratios. 



The ^'^Si(-'He, ny) S Experiment - Phase II Theory 



Ey (keV) 


JM 


* Jm 






Mult." 


5 


Rdco 


A/ 

ijm^ if) 


Mult." 

Om^ jf) 


AJ 

(ji^ jm) 


Mult." 

iji^ jm) 


Rdco 


1194.0(1) 


2- 


->2- 


^0 


0.5* 


M\/E2 


0.16* 


0.92(4) 


2 


E2 





D 


1.0 


1456.5(3) 


0- 


-^2- 


^0 


0.3^^ 


E2 


Qd 


0.94(9) 


2 


E2 


2 


Q 


1.0 


1466.2(3) 


1 - 


-^2- 


^0 


0.3^^ 


Ml, £2 


-0.09(3)'' 


0.40(8) 


2 


E2 


1 


D 


0.5 


2477.3(3) 


3- 


-^2- 


^0 


0.3'' 


Ml, £2 


0.73(9)'' 


0.37(4) 


2 


E2 


1 


D 


0.5 


2921.4(4) 


4- 


-^2- 


^0 


0.3^^ 


E2 


Qd 


0.99(11) 


2 


E2 


2 


Q 


1.0 



"Transition multipolarity, D and Q refer to dipole and quadrupole, respectively. 
*Determined experimentally from angular distribution measurements. 

■^The alignment factor of a/ ; = 0.3 is usually adopted when no experimental information is available for this parameter. Since y-ray angular distribution 
measurements were only obtained for the 2210.6-k:eV and 1 194-keV y-rays, we have assigned oj j = 0.3 for all other y-rays of ^^"S. 
''From selection rules. 

''This mixing ratio was adopted from the mirror transition (see Ref. [40]). 



it can be used as a guide for determination of the spins of the 
associated states involved in the transitions. 

For a jf j,n — > ji cascade, where jf, j,„ and ji are the 
spins of the final, intermediate and initial states, respectively, 
the DCO ratio is generally defined as [55]: 

^e'(Gate^') 
Rdco — — yTt ' 

where 9i and 02 are the angles with respect to the beam axis 
at which detectors 1 and 2 are placed, respectively; / is the in- 
tensity; and Yi and 72 are transitions observed in coincidence, 
which originate from the j„, — > y, and jf — > y„, decays, respec- 
tively. Yi transition is the one on which the coincidence gate 
is placed. 

The theoretical DCO ratios are given in Table V for the 
cases where the Y2 transition from the jf — y„, decay is a 
pure transition. If on the other hand, Y2 is a mixed transition, 
the theoretical DCO ratio is expected [54] to differ from what 
is listed in Table V. The significance of such a difference de- 
pends upon the severity of the dipole-plus-quadrupole admix- 
ture of the Y2 transition [44] . 

For our y-Y angular correlation measurement, detectors 1 
and 2 were positioned, with respect to the beam axis, at 90° 
and 135°, respectively, and on opposite sides with respect to 
the beam line. The — > 0,*" transition in ^"S was measured 
with detector 1 and after gating on this transition, the higher 
lying transitions were observed in the coincidence spectrum 
measured by detector 2. Both detectors were placed as close 
to the target as possible, i.e., 3 cm and 7 cm away from the 
target, respectively. 

The y-Y angular correlations of ^"S y-rays were determined 
by measuring the DCO ratios for each ^"S y-ray that was ob- 
served at both angles. Since the statistics under the 1283.4- 
keV, 1405.1-keV, and 2599-keV y-ray transitions correspond- 
ing to the 3| 2j, 2| 2j and 2j -> 2| decays in ^"S, 
respectively, are too poor, the DCO ratio could not be deter- 
mined for these transitions. 

The experimental DCO ratios of all other transitions of 



were determined after a gate was set around the stretched 
quadrupole transition with 2210.6-keV energy (2j^ 0^), ob- 
served in the singles y-ray spectra at both angles, to obtain 
the corresponding coincidence spectra. The peaks of interest 
in the coincidence spectra were then fitted, and their yields, 
corrected for detector efficiencies, were obtained and used to 
calculate the DCO ratios via Eq. (9). The results are given in 
Table V. 

The y-y directional correlations of y-ray transitions in ^"S, 
and the spin-parity assignments of ^"S states with E^^ < 5. 136 
MeV from our (p,t) measurements are discussed below. 

(2.4). Spin-parity assignments 

Prior to discussing the spin-parity assignments, it should 
be noted that in the following discussion, the energies of the 
adopted y-ray transitions are corrected for '^"S recoil energies 
(see Table III and Fig. 6). 

• The 2210.7-keV y-ray Transition: corresponds to the 
2210.7-keV ^ ground state decay transition. The 2210.7- 
keV state was observed only at 62° during phase I of the 
^^S(/9,r)^*'S experiment, thus no J'^ assignment is available 
from that experiment. However, our present y-ray angular dis- 
tribution parameters for this transition (see Table IV) confirm 
r = 2+. 

• The 1194-keV y-ray Transition: corresponds to the 
3404.7-keV 2210.7-keV decay transition. Due to the lack 
of triton angular distribution data from our ^^S(/7,f)^"S exper- 
iments for the 3404.7-keV state, no conclusive spin-parity as- 
signment was obtained for this state from those experiments. 
Nevertheless, the 7" assignment for the 3404.7-keV state is al- 
ready established as 2+ from various previous measurements, 
e.g., Ref. [13], and the results of our y-ray angular distribution 
measurements for the 1 194-keV y-ray transition agree with a 
A/ = transition from a 7" = 2+ state (see Table IV). IVIore- 
over, our experimental Rdco ratio for the 2^ ^ Qi\ 
cascade agrees with the theoretical ratio within 2a, and is con- 
sistent with an Ml transition with a small E2 admixture for the 
1 194-keV y-ray transition for the 2 J — > 2'^ decay. 
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• The 1456.5-keV y-ray Transition: corresponds to the 
3667.2-keV 2210.7-keV decay transition. The 3667.2-keV 
state could not be resolved in our ^^Sip,t)^^S experiments, 
and thus no information on its energy or spin-parity is avail- 
able from those experiments. The present experimental and 
theoretical Rdco ratios for the — > 2j'^ — > Oj^ cascade are 
consistent with unity, suggesting that the transition from the 
3667.2-keV state to the 2210.7-keV state has the same mul- 
tipolarity as that of the decay of the 2210.7-keV state to the 
ground state (see Ref. [55] and Table V). This implies that the 
1456.5-keV y-ray is a pure quadrupole transition. Therefore, 
we confirm the assignment of = 0+ for the 3667.2-keV 
state, because from the mirror nucleus no other possibilities 
are expected in this energy range for a A/ = 2 transition corre- 
sponding to the 3667.2-keV — > 2210.7-keV decay transition. 

• The 1466.2-keV y-ray Transition: corresponds to the 
3676.9-keV — > 2210.7-keV decay transition. From our 
■^^S{p,t)^^S experiment, we obtained an energy of 3681(3) 
keV, consistent with the 3676.9(4) keV obtained from our in- 
beam y-ray spectroscopy experiment within 2a. The present 
triton angular distribution data for the 3681-keV state agree 
with both /"^ = 0+ and 7" = 1+ (see Fig. 3). Previous mea- 
surements [13, 16] have assigned a ^ = 1+ to this state. Ac- 
cording to Table V, the theoretical Rdco is expected to be 
0.5 if the 1466.2-keV y-ray transition is a stretched dipole 
(£1 or Ml transition with 5 = 0) A/ = 1 transition from a 
state with = 1+ or = 3+. If, on the other hand, the 
aforementioned transition is a mixed dipole-plus-quadrupole 
instead of a stretched dipole, the theoretical Rdco should dif- 
fer from 0.5 [54]. Considering the 7'' — > 2+ — > 0^ cascade 
as the 3676.9-keV — > 2210.7-keV ground state decay tran- 
sitions, our previous discussion implies that the 3676.9-keV 
state could either be the 1 or 3| state in ''^S. A J'^ = 0+,l + 
doublet is thought [15, 16] to exist in 3.6 - 3.8 MeV 
region in ^"S. Being very close in energy to the 3667.2-keV 
state, the 3676.9-keV state must be the 1+ member of the 
aforementioned doublet, now that we have confirmed the for- 
mer as the 0+ member Our experimental Rdco ratio for the 
3676.9 keV — > 2210.7 keV — > ground state cascade is slightly 
lower than 0.5 (see Table V), which implies that the 3676.9- 
keV state is most likely the l| state of ^"S and the 3676.9 
keV — >■ 2210.7 keV decay transition is a likely an Ml tran- 
sition with a small E2 admixture. We could not determine 
the mixing ratios of any of the transitions observed via the 
present y-y directional correlation measurements. Therefore, 
we have adopted the mixing ratio of -0.09(3) [40] (from the 
mirror transition) for the 1466.2-keV y-ray transition of ^''S. 
In conclusion, we suggest a 7" = 1+ for the 3676.9-keV state 

of 30s. 

• The 2477. 3 -keV y-ray Transition: corresponds to the 
4688.1-keV — > 2210.7-keV decay transition. The present tri- 
ton angular distribution for the 4688.1-keV state is consistent 
with a = 3+ assignment (see Fig. 3). Moreover, the decay 
branches of the 4688-keV state, observed in our in-beam y-ray 
spectroscopy experiments, also agree with those of the mirror 
state [24] assuming that the 4688-keV state is the 3 state of 
^"8. Our experimental DCO ratio for the 3| ^> 2[^ — > 0,^ cas- 
cade is significantly different from the theoretical Rdco — 0-5 



(see Table V). Therefore, based on the previous discussion, 
we expect the 2477.3-keV y-ray to be a A7 = 1 mixed MIIE2 
transition from a = 3+ or a = 1+ state. According to the 
mirror states in ^^Si [40], only one 7" = 1 + state is expected in 
this energy range, and that is most likely the 3676.9-keV state. 
These arguments suggest that the 4688.1-keV state is the 3]^ 
state of Therefore, our experimental /?£)co ratio also sup- 
plements the other present results with regards to the value 
of the 4688-keV state. We have adopted the mixing ratio of 
the mirror transition (5 = 0.73(9) [40]) for the 2477.3-keV y- 
ray due to the lack of knowledge of its own mixing ratio. We 
conclude that the 4688-keV state is the mirror to the 2>'\ state 
in ^^Si at 4831 keV [40]. Thus, the 4688-keV level in ^"8 is 
the 3j" astrophysically important state predicted by lliadis et 
al. [9]. 

• The 2599.1 -keV y-ray Transition: corresponds to the 
4809.8-keV 2210.7-keV decay transition, which is a 
very weak transition observed at 135° in the present 
^^8i(3He,ny)30S experiment. Therefore, no experimental 
Rdco ratio could be obtained for this transition. The present 
triton angular distribution data agree with both 7" ~ 2+ and 
3+ (see Fig. 3) but the former is a better fit. The y-ray branch- 
ing ratios for the y-decay of the 4809.8-keV state to the 2'^ and 
2 J states in were measured at 90° [24], and were in good 
agreement within their uncertainties with those of the decay of 
the 4810-keV state in ^^Si to its 2'^ and 2t lower-lying states 
(also see § 11 B 2). Also, in a recent shell model calculation for 
the sd-shell in A = 30 nuclei using the U8D Hamiltonian with 
inclusion of a charged-dependent term [56, 57], the energy of 
the 2'^ state in was derived to be near 4800 keV, while 
that of the 3 state was calculated to be near 4700 keV. These 
results altogether strongly support a 7'^ = 2+ assignment for 
the 4809.8-keV state (mirror to the 2| state at 4810-keV in 
3^81 [40]), making it the next astrophysically important state 
predicted by lliadis et al. [9]. 

• The 2921.4-keV y-ray Transition: corresponds to the 
5132.3-keV 2210.7-keV decay transition. In the shell- 
model analysis by Wiescher and Gorres [11], they concluded 
that there are most likely at least two levels with energy near 
5 MeV; a 4+ near 5.1 MeV and a 0+ near 5.2 MeV. Kuhlmann 
et al. [16] observed a state at 5136(2) keV, and concluded that 
this level is most likely a 4+ state. In Ref. [13], a state was 
observed at 5168(6) keV. The triton angular distribution data 
in that work could not be fitted with a single angular momen- 
tum transfer, which suggested that the latter state was an un- 
resolved doublet consisting of a 4+ and a 0^ state. Our triton 
angular distribution data are best fitted with / = 4 transfer 
Our experimental Rdco ratio for the 4'^ ^ 2^ ^ 0^ cascade 
is consistent with the theoretical ratio given in Table V under 
the assumption that the 5132.3-keV — > 2210.7-keV transition 
is a stretched quadrupole with A7 = 2. This indicates that the 
5132.3-keV state, observed in our in-beam y-ray spectroscopy 
experiments, is either the 4j^ or the 0^ state of ^"8. The for- 
mer is much more probable because a comparison with the 
mirror transitions in ^^Si reveals that the 4j^ level at 5279.37 
keV in ^"81 decays with a 100% branch to the first excited 2'^ 
state [40]. This is consistent with what we observe for the 
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Table VI: Energy levels of ■'^S from this work with 
^^S{p,t)^^S measurements are not shown here. 



< 6 MeV . The energies of the states used as internal calibration energies in our 



Present Work 



Present Work 



Adopted Energy 
(keV) 



(keV) 



Ex (keV) 




Ex (keV) 




Ex (keV) 






2208(3) 




g.s. 
2210.7(3) 


2+ 


g.s. 
2210.6(3) 


^TX 

0+ 
2+ 








3404.7(3) 


2+ 


3403.6(6) 


2+ 








3667.2(4) 


0+ 


3667.0(5) 


0+ 




3681(3) 


(l+,0+) 


3676.9(4) 


1+ 


3677.0(4) 


1+ 




4688(2) 


3+ 


4688.1(4) 


3+ 


4688.1(4) 


3+ 


293.2(8) 


4812(2) 


2+ 


4809.8(5) 




4809.8(6) 


2+ 


414.9(9) 




(4+)* 


5132.3(5) 


(4+) 


5132.6(8) 


(4+) 


737.7(11) 


5225(2) 


(0+) 






5221(2) 


(0+) 


826(2) 


5315(2) 


(3-, 2+) 






5314(4)^^ 


(3-) 


919(4) 


5393(2) 


3+ 






5391(2) 


3+ 


996(2) 


5849(2) 


(l",2+,4+) 






5847(2) 


(2+) 


1452(2) 


5947(2) 








5946(3) 


(4+) 


1551(3) 



"Er = £.1 — Q, where £",- is the resonance energy, ZJj is the weighted average excitation energy and Q is the proton threshold of the P(p,y) S reaction (4394.9 
keV). Those excitation energies for which no resonance energy is reported correspond to the bound states of ^"S. 
''The con'esponding energy (5 1 36(2) keV) was used as internal calibration energy, and is thus not reported here. 
''This state is most likely the 5288-keV state observed by Yokota et al. [17], which was assigned to be the 3J" state in '^"S. 



5132.3-keV state in ^"S, as well as what was observed for the 
same state in Ref. [16]. If the 5132.3-keV state were the O3" 
state, based on its decay scheme in the mirror nucleus, we 
would have expected to observe other decay branches from 
this state with comparable strengths, in addition to the 292 1 .4- 
keV y-ray transition [40]. From these arguments, we tenta- 
tively assign 7" = 4+ to the 5132.7-keV level of ^"S. 

In the following subsection, the spin-parity assignments for 
a few other ^"S states with ^ 6 MeV are discussed. 



(2.5). Adopted energy levels in ^^S 

Table VI presents the combined results of both phases of 
both our experiments on ^"S excitation energies below 6 MeV, 
and the corresponding recommended spin-parity assignments. 
The adopted energies in Table VI are the ^"S weighted average 
excitation energies over all independent measurements in the 
literature, including the present work. States used as internal 
calibration energies were excluded in the calculations of the 
adopted energies. In a few cases where the uncertainty in the 
weighted average was smaller than the smallest uncertainty 
in the measured excitation energies, the latter was adopted as 
the final uncertainty only if the energy was measured in fewer 
than 4 independent measurements [44]. 

From our {p,t) measurements, a unique spin-parity assign- 
ment could not be determined for the 1452-keVand 1551-keV 
resonances. To calculate the ^^P(/7,y)^''S reaction rate, a 
value had to be assumed for each of these resonances. 

A tentative value of 1^ was assigned [13] to the 1452- 
keV resonance (see Table VI), but due to poor statistics for this 
particular resonance, / = 2 or 3 transfers were not excluded. 
In our {p,t) measurements, the 7" value for the 996-keV res- 
onance fits best with a 3+ assignment, and we have assigned 



the 91 9-keV resonance to be the 3 J state (see § II A 2). Hence, 
we have tentatively assigned the 1452-keV resonance to be the 
2| state in ^"8. 

The 1551-keV resonance has only been observed in the 
measurement of Ref. [18] and in our ^^S(/9,f)^''S measure- 
ments. However, the data obtained in these measurements 
were not enough to assign a conclusive value to this res- 
onance. From the results of a recent shell model calcula- 
tion [56, 57], the energies of the ^2 states in ^"S and its mirror 
nucleus - ^"Si - are almost identical to each other. The exci- 
tation energy of the 4j state in ^^Si is 5950.73(15) keV [40]. 
On the other hand, the weighted average energy between the 
results of Ref. [18] and our ^^S{p ,t)^^S measurements for the 
corresponding state in ^"8 is 5946(3) keV. Therefore, we con- 
cluded that this latter state is most likely the mirror to the A2 
state in ^^Si at 5950.73(15) keV. 

The states presented in Table VI are the only ones that could 
play a crucial role in determination of the ^^P(/7,y)^''S reac- 
tion rate in the temperature range characteristic of explosive 
hydrogen burning (0.1 GK <T < 1.3 GK). The excited states 
whose energies are below 4.5 MeV become important in de- 
termining the non-resonant contributions to the '^P(p,Y)^''S 
reaction rate. 



III. THE -''P(p, Y)^''S REACTION RATE 

To obtain the non-resonant contribution to the ^^P(/?,y)^''S 
reaction rate, one has to determine the astrophysical S-factor, 
S{E), from: 

S{E) w 5(0) + S'{0)E + ^S"{0)E^, (10) 
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where the primes indicate derivatives with respect to E. S{E) 
can be integrated to give the non-resonant reaction rate: 

(11) 

where Na < OD > is the reaction rate, Na is Avogadro's num- 
ber, /J is the reduced mass, k is Boltzmann's constant, T is the 
temperature, E is the center-of-mass energy in keV, and Eg is 
the Gamow energy (in MeV). 

The '^P(p,y)^''S direct capture (DC) reaction rate to all 
bound states, including the ground state, was calculated as- 
suming proton transfer into 2s and Id final orbitals. For each 
final state, the S-factor was calculated by taking into account 
the £1 and Ml nature of the transitions, which were then 
weighted by the corresponding spectroscopic factors deter- 
mined from those of the mirror states [58]. The weighted S- 
factor contributions from each state of ^"S were then summed 
to derive the total S-factor as a function of proton bombarding 
energy for each transition multipolarity. 

The S-factor was then fitted with a polynomial of the form 
given in Eq. (10) to determine the fit parameters, i.e., 5(0), 
S'{0) and 5"(0). An uncertainty of 40% for the direct capture 
S-factor is adopted following the approach of Ref. [59]. 

With increasing center-of-mass energy, resonances become 
important, and therefore the non-resonant S-factor in Eq. (10) 
is truncated at the so-called cutoff energy, after which the di- 
rect capture S-factor deviates from the total astrophysical S- 
factor. The cutoff energy was chosen [59] at ^1000 keV for 
the ^^P(/7,y)'"'S reaction rate. 

To calculate the resonant contributions to the rate, the pro- 
ton widths were determined using the expression: 

rp = 2^PiC^SQ%, (12) 

where /j is the reduced mass, P/ is the barrier penetrability 
(calculated using ro = 1.25 fm) for orbital angular momen- 
tum I, a = ro(Ay^ +aI/^) is the interaction radius in terms of 
target and projectile mass numbers (A, and Ap, respectively), 
C and S are the isospin Clebsch-Gordan coefficient and spec- 
troscopic factor, respectively, and 9jp is the observed dimen- 
sionless single-particle reduced width. 

The 9jp factors were estimated using Eq. (11) together with 
Table 1 of Ref. [60]. The only exceptions were the 4+ states 
corresponding to the resonances at 737.7 keV and 1551 keV. 
The reduced widths of these resonances could not be deter- 
mined from the approach of Ref. [60], which is limited to 
single-particle states in the sd - fp shells. Consequently, 
9?_ < 1 is assumed for these states. 

Spectroscopic factors were determined from neutron 
spectroscopic factors of the mirror states measured with 
^'^Si{d,p)^°Si [58]. The mirror levels in 30si corresponding 
to the resonances of ^°S at 737.7 keV, 826 keV and 1551 keV 
were populated very weakly in the measurement of Ref. [58], 
and thus no C^S values could be determined experimentally 
for these levels. Hence, an upper limit of C^S < 0.01 is 
adopted for these states, based on the sensitivity for the ex- 



traction of small spectroscopic factors. Following the proce- 
dure of Ref. [59], the uncertainties in the proton widths were 
estimated to be 40%. 

To determine the y-ray partial widths (Fy), the correspond- 
ing widths of the mirror states in ^''Si were calculated from 
measured half-lives, branching ratios, multipolarities, and 
mixing ratios [40]. For the cases where mixing ratios of the 
transitions of interest in 30si have not been determined exper- 
imentally or theoretically, we have assumed that such tran- 
sitions are pure, with multipolarities assumed to be the domi- 
nant multipolarity of the actual mixed transition. These widths 
were then scaled to account for the energy difference between 
each mirror pair, assuming similar decay branches and re- 
duced transition probabilities. 

Only an upper limit is known for the half-life of the 2^ state 
in 30si. Hence, the aforementioned method of calculation of 
the y-ray partial width is not possible for the corresponding 
mirror state in ^''S. Therefore, we considered the 0.012 eV 
value from Ref. [59] for the total y-ray width of the 2^ res- 
onance in ^'^S, and scaled it to account for the differences in 
the measured energies. Following the procedure discussed in 
Ref. [59], the uncertainties in y-ray widths are assumed to be 
50%. 

The ^^P(p,y)^''S reaction rate was calculated using the 
Monte Carlo method presented in Refs. [12, 59, 61] using 10"* 
random samples. The input file required for calculation of 
the ^^P(/7,y)^*^S Monte Carlo reaction rate, including the reso- 
nance parameters, is given in Ref. [25]. The numerical values 
of the 2^P(/7,y)^"S rate ai-e given in Table Vll. 

Since the experimental input is truncated at 6 MeV (see 
Table VI), there exists a temperature denoted by Thatch [12], 
above which the nuclear physics input is insufficient to calcu- 
late a reliable reaction rate. For our case, r„,„,c/j was found to 
be 3 GK. The reaction rate above this temperature was there- 
fore obtained as follows: the theoretical ^^P(/7,y)^"S rate as a 
function of temperature was calculated based on the Hauser- 
Feshbach statistical model [62] using the non-SMOKEr'^^^ 
code [63]. Then, these rates were scaled and normalized to 
the experimental Monte Carlo rate at Thatch = 3 GK. 

Fig. 1 1 compares the contributions of the direct capture rate 
(DC Rate) and those of the resonances listed in Table VI to the 
total ^^P(/5,y)^''S thermonuclear reaction rate. The direct cap- 
ture rate dominates the total rate for T < 0.08 GK, whereas at 
higher temperatures characteristic of explosive nucleosynthe- 
sis in novae, the total reaction rate is dominated by a single 
3j^ resonance at 293.2 keV in the range of 0.09 - 0.3 GK. 
The 414.9-keV resonance with 7" = 2^^ is the main contrib- 
utor to the total rate from 0.35 GK to 2 GK. The 996-keV 
resonance becomes important at temperatures higher than 2 
GK, which are beyond the temperature range of interest to this 
work. The other resonances, including those for which only 
an upper limit proton partial width is known, do not contribute 
significantly to the "^P(/:',y)^*'S total rate in the temperature 
range of interest. 

Fig. 12 compares our Monte Carlo rate for the ^^P(p,y)^''S 
reaction with that of Ref. [12], where the energies of the 3j^ 
and 2^ states of ^"S were assumed to be 4704(5) keV [13] and 
4888(40) keV [9], respectively. 
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Table VII: Total Monte Carlo rate for the P(p, y) S thermonuclear reaction. See text for details. 



T(GK) Low Rate Median Rate High Rate T (GK) Low Rate Median Rate High Rate 



0.010 


5.06x10- 
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important resonances (corresponding to the 3) and 23 states 
in ^"S) were determined theoretically based on the IMME, 
since none of these resonances were observed at the time. 

For our present rate, the energy of the resonance corre- 
sponding to the 2^ state of ^'^S is determined experimentally. 
Thus, its uncertainty of 0.9 keV is reduced by a factor of ~44 
with respect to the theoretical estimate of 40 keV adopted in 
Ref. [12]. Furthermore, the 2-keV uncertainty in the energy 
of the resonance corresponding to the 3^ state in ^"S, mea- 
sured in this work, is also reduced by 40% with respect to the 
5 keV measured in Ref. [13] that is used to derive the rate in 
Ref. [12]. 

Therefore, since these two resonances together dominate 
the total rate over 0.08 < T < 2 GK, the reductions in their 
associated uncertainties also reduce the uncertainty in the total 
reaction rate (see the bottom panel in Fig. 12). For example, 
at r = 0. 1 GK, where the uncertainty in both our rate and that 
of Ref. [12] is maximum, the Na < OV >hi,JNA < OV >i„„ ra- 
tio from our Monte Carlo rate is 72% smaller than that of the 
Monte Carlo rate reported in Ref. [12]. 



IV. NOVA ISOTOPIC ABUNDANCES 

In § I, it was emphasized that ^^P(p,Y)^*'S is one of the two 
reactions that are thought to affect the silicon isotopic ratios 
in nova ejecta. 

To investigate the impact of the updated ^^P(/7,y)^*'S rate 
on the isotopic abundances of silicon synthesized in classical 
novae, we have computed three different models of nova out- 
bursts, with identical input physics except for the adopted 
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Figure 11: (Color online) Resonant and DC contributions to the 
^^P(P,Y)^*'S Monte Carlo rate as a function of temperature. Abbre- 
viations are as follows: DC: Direct Capture; UL: Upper Limit. The 
latter is used for those resonances for which the proton partial width 
is estimated to be an upper limit. 



Both rates shown in Fig. 12 are calculated using the Monte 
Carlo technique. Our lower resonance energies, relative to 
those adopted in Ref. [12], for the two astrophysically im- 
portant resonances, cause our median rate to be up to 2.3 
times larger (at Tg =0.1) (see top panel in Fig. 12) and up to 
11.5 times larger (at Tg = 0.1) than the recommended rate of 
Ref. [9]. For this last rate, the energies of both astrophysically 
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Figure 12: (top) The ratio of our Monte Carlo low (lower thin line), 
median (middle thick line) and high (upper thin line) rates to the 
Monte Carlo median rate of Ref. [12]. Our median rate is 2.3 times 
larger than that of Ref. [12] at T = 0.1 GK. (bottom) The uncer- 
tainty bands corresponding to A'^ < OD >i,igi,/A'^ < OV >mLdiiin and 
Na < OV >]„JNa < OV >mcdiim from our Monte Carlo rate (solid lines) 
compared to those of Ref. [12] (dashed lines). At T = 0. 1 GK, the 
ratio of the A'^ < av >i,i„i,/A'^ < OV from our Monte Carlo rate is 
72% smaller than that of the Monte Carlo rate reported in Ref. [12]. 



P(/7,Y) S rate. Results from our nova nucleosynthesis sim- 
ulations are presented next. 



A. Nova simulations 

Three nova nucleosynthesis models were computed with 
the Lagrangian one-dimensional (spherically symmetric) full 
hydrodynamic and implicit code called SHIVA. Detailed in- 
formation about this code is provided in Refs. [64, 65]. 

SHIVA simulates the evolution of nova outbursts from the 
onset of accretion to the explosion and ejection of the nova 
ejecta. The hydrodynamic code is coupled directly to the nu- 
clear reaction network. Thus, in the present work, the explo- 
sion simulations at each stage are complemented with detailed 
nova nucleosynthesis calculations using the most updated re- 
action rate libraries. 



As pointed out in Ref. [64], the material is dredged up on 
short timescales from the outermost shells of the CO- or ONe- 
rich core to the surface of the white dwarf by convective mix- 
ing processes. Nuclear reactions in stellar environments are 
sensitive to the temperature, and thus the ejected abundances 
of fragile nuclei that would have been destroyed if they had 
not been carried to higher and cooler layers, are increased 
by considering the convection process during the evolution 
of the nova outburst. This, in turn, makes the present sim- 
ulations more realistic and suitable for defining absolute iso- 
topic abundances resulting from nova nucleosynthesis than the 
previous post-processing nucleosynthesis simulations used in 
Refs. [1, 13], where the nucleosynthesis is decoupled from the 
hydrodynamics of the outburst. 

The absolute abundances observed in nova ejecta or in 
presolar grains of potential nova origin provide strong con- 
straints for improvement of nova simulations. Thus, a more 
precise set of constraints can be obtained if predictions on spe- 
cific isotopic abundances are available. 

For the present full hydrodynamic simulations, the ther- 
modynamic profiles are identical to those of hydrodynamical 
simulations, given in Ref. [64], for a massive ONe nova with 
a 1.35 Mq underlying white dwarf. Such an extreme white 
dwarf is adopted because a CO white dwarf shows limited ac- 
tivity in the Si-Ca mass region. This, in turn, is due to very 
little, if any, Ne, Mg and Si seed nuclei available in the outer 
core of a CO white dwarf, and the lower temperature achieved 
in a CO nova outburst [7]. Thus the nucleosynthesis of silicon 
isotopes in CO novae, with even the most massive underlying 
white dwarf, is negligible. 

An accretion at a rate of = 2 x lO^'^MQ/yr of solar- 
like matter onto a 1 .35 Mq ONe white dwarf is assumed in all 
three present models. 50% enrichment by the white dwarf's 
core material is adopted for the accreted matter to mimic the 
unknown mechanism responsible for the enhancement in met- 
als, which ultimately powers the explosion through hydrogen 
burning [7]. The initial abundances of the seed isotopes used 
in the present simulations are given in Ref. [7]. The impact of 
the new solar metallicity [66] (decreased by about a factor of 
2) on the overall results presented here has been tested and is 
insignificant. 

In addition to hydrodynamics, a reaction rate network in- 
cluding 370 nuclear reactions involving 117 isotopes rang- 
ing from 'H to '^'^Ti is used. Monte Carlo reaction rates are 
adopted from the most updated compilation of Ref. [12] with 
additional reactions selected from the reaction rate library of 
Iliadis (2005). The only exception is the ^'^P(p,y)^^S reaction, 
whose rate is chosen (one at a time for each of the three mod- 
els) from the present work, as well as from Refs. [9, 12] for 
comparison. 

These rates are corrected for the stellar enhancement factors 
to allow for the increase in reaction rates associated with par- 
ticipation of excited states of nuclei in the reactions. Lastly, 
the impact of the ^^P(/9, y)^''S stellar reaction rate on nova nu- 
cleosynthesis was compared for the three different reported 
rates: the recommended classical rate from Ref. [9], here- 
after model A; the median Monte Carlo rate from Ref. [12], 
hereafter model B; and the high Monte Carlo rate from this 
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work, henceforth model C. The main distinctions in the three 
2^P(/7,Y)^*^S rates used in the present nova simulations arise 
from different input energies and uncertainties for two reso- 
nances corresponding to the 3 and 2'^ states of ^''S (see § 111 
for discussion). 

The selection of the high Monte Carlo rate from this work 
instead of the median rate is to account for the highest possible 
effect of the new rate on the abundances of elements synthe- 
sized in novae. While our median rate is 2.3 times larger (at 
0.1 GK) than that of Ref. [12] (see top panel in Fig. 12), the 
present high rate is a factor of 3.5 and 17 larger (at 0.1 GK) 
than the median rate of Ref. [12] and the recommended rate 
of Ref. [9], respectively. 



B. Results 

To assign different weights to individual shells of the under- 
lying white dwarf, the isotopic abundances obtained from the 
three aforementioned hydrodynamic nova simulations were 
averaged over mass within each shell. The total ejected en- 
velope mass for each of the three models is 4.55 x 10^^ 
Mq. The resulting mean abundances (in mass fractions) in 
the ejected envelope shells for models A to C are given in Ta- 
ble Vlll for a selection of the stable isotopes in the Si-Ca mass 
region, whose abundances (in mass fractions) are greater than 
or equal to 10^^. Those stable isotopes not included in Ta- 
bles Vlll did not change significantly between models. 

For the stable isotopes with 14 < Z < 20 which are products 
of the decays of the short-lived radioactive species, a compar- 
ison was made between the mean abundances obtained from 
model C and those obtained from models A and B. With re- 
spect to models A and B, the largest abundance change ob- 
served from the results of model C is a 6% decrease in the 
abundance of ^^Si. This percentage difference is defined to 
be: [(new value — old value) ^ old value], where the "new" 
value is an isotopic abundance or ratio resulting from model 
C, and the "old" values are those resulting from models A or 
B, whichever gives a higher percentage difference. A nega- 
tive (positive) percentage difference indicates that the isotopic 
abundance or ratio resulting from model C is decreased (in- 
creased) with respect to that obtained from model A or B. 

Therefore, changing the ^'^P(p,j)^'^S rate seems to have 
only a small effect on the abundances of isotopes with A « 30 
produced in a nova outburst. However, because of the reduced 
uncertainty in the updated ^'^P{p,"f)^'^S rate, we are now more 
confident in the reliability of the isotopic abundances obtained 
using model C. 

The abundance of each stable isotope alone does not pro- 
vide much useful information. Instead, to compare the iso- 
topic abundances obtained from nova simulations with those 
observed in presolar grains, one has to investigate an isotopic 
abundance ratio. For example, the silicon isotopic ratios mea- 
sured in presolar grains are usually expressed as [7]: 
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Table Vlll: Selected mean composition of nova ejecta (in mass frac- 
tions, for the Si-Ca isotopes) from models of nova explosions on 1.35 
Mq ONe white dwarfs. The only difference between models A, B 
and C is the 25p(/7,Y)3"S rate used. 



Isotope 



Hydrodynamic Model 





A 


B 


C 




Ref. [9] 


Ref. [12] 


present work 




(recommended) 


(median) 


(hieh) 


28Si 


3.08xlO-°2 


3.08xlO-°2 


3.08xl0-"2 


29si 


2.38xlO"°3 


2.39 xlO"*" 


2.24x10"'" 


3"Si 


1.54xl0-°2 


1.54xl0-*'2 


1.51xlO-"2 


31p 


8.71xlO"°3 


8.73x10"*" 


8.61 xlO"'" 


32s 


5.27xl0"°2 


5.27xl0"°2 


5.30xlO""2 


"s 


8.02x10"""* 


8.01 X lO"*'"* 


8.17x10""'* 


34g 


3.63 xlO"""* 


3.63xl0-°4 


3.71x10""'* 


35C1 


3.85xlO-°4 


3.85xlO-°4 


3.95xlO-"4 


36^ 


5.14x10"°-** 


5.14x10"°-'* 


5.29x10""^ 


3**Ar 


2.19x10""^ 


2.19x10"°^ 


2.21 xlO"'" 



Table IX: Deviations (in permit) from solar abundances in simulated 
and measured nova silicon isotopic abundances. Models A to C are 
explained in the text, and are obtained from hydrodynamic simula- 
tions of classical nova outbursts. The measured values (the first four 
rows) are for SiC presolar grains reported in Refs. [6, 7]. 



Grain 


5(2'^Si/2«Si) 


5(3"Si/2'^Si) 


Hydrodynamic 




(%c) 


(%o) 


Model 


AF15bB-429-3 


28 ±30 


1118 ±44 




AF15bC- 126-3 


-105 ± 17 


237 ± 20 




KJGM4C- 100-3 


55 ±5 


119 ± 6 




KJGM4C-311-6 


-4 ±5 


149 ±6 






527.1 


13970 


A 




533.5 


13970 


B 




437.3 


13678 


C 



where 8 represents deviations from solar abundances in per- 
mil, and the adopted numerical values for the solar silicon 
isotopic ratios are [67] (p. 130): (2''Si/2^Si)Q = 0.0506 and 

(3"Si/2^Si)Q = 0.0334. The deviations from solar abun- 
dances are computed for silicon isotopic abundance ratios ob- 
tained from models A, B and C, and the results are shown 
in Table IX, along with ^'^ 30|§jy28^j ratios measured [6] from 
some SiC presolar grains with proposed classical nova pater- 
nity. 

As seen in Table IX, the theoretically predicted 8-values 
are much larger than the measured counterparts. Overall, 
however, regardless of the ^^P(p,y)^^S reaction rate used, the 
2^Si/2^Si ratio in the ejecta resulting from the simulations 
is only slightly higher (~1.5 times larger, see Table Vlll) 
than the solar value. Using the measured 5(2^Si/2^Si) values, 
given in Table IX, as inputs to Eq. (13), we extract a mea- 
sured 2^Si/2^Si ratio that varies between a factor of 0.9 - 1.1 
times the solar ratio, and thus is again only slightly lower or 
higher than the solar value. Therefore, even though the new 
2^P(;5,y)^"S rate does not significantly improve the theoretical 
5-values, the simulated signatures are qualitatively consistent 
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with the ^Si/ Si ratios measured in presolar grains identi- 
fied to have a nova origin. In other words, the simulated and 
measured 5-values both show enhancements in the same di- 
rection. 

On the other hand, the ^''Si/^^Si ratio in the ejecta result- 
ing from the simulations is much higher (~15 times larger) 
than the solar value (see Table VIII), such that the classical 
nova ejecta resulting from the hydrodynamic models is sig- 
nificantly enriched in 30si. The simulated and measured val- 
ues again are in qualitative agreement with each other, i.e., 
enhanced in the same direction, but the magnitudes of the en- 
hancements are not in agreement. 

Our results support the indication that in order for the 
models to predict the ^"Si/^^Si ejecta ratio that quantitatively 
matches the grain data, one has to assume a mixing process 
between material newly synthesized in the nova outburst and 
more than 10 times as much unprocessed, isotopically close 
to solar, material before the process of grain formation [5, 7]. 
The details of the ejecta dilution and the grain formation pro- 
cesses are still unknown. 

In addition to invoking the mixing with solar composi- 
tion material, an increase in the ^'^P{p,y)^^S reaction rate also 
helps reduce the ^"Si/^^Si ratio by moving the nucleosynthesis 
flow away from ^"P toward the heavier isotopes. A decrease 
in the abundance of ^"P consequently reduces that of ^''Si pro- 
duced from 3''P(p+)3''Si. The rate of the ^^P(p,y)^^S reaction 
has been evaluated in Refs. [8, 68, 69] and more recently in 
Ref. [70]. This last rate is found to be ^10 times greater, at 
T 0.25 GK, than the lower limit set in Ref. [69]. A factor 
of ~^10 increase in the ^''P(p,y)-''S reaction rate results in a 
typical factor of -^10 reduction in the expected abundance of 
^"Si [70]. This new information may now help to better con- 
strain the dilution process in new nova model predictions. 

In comparison with the high Monte Carlo rate from the 
present work, the present median and low Monte Carlo rates 
show smaller deviations with respect to the median rates of 
Refs. [9, 12]. Therefore, we did not extend our investigation 
to study the effects of these rates on the nova yields. 

V. CONCLUSIONS 

The ^^P(p,y)^^S reaction rate at the temperature range of 
0. 1 < r < 1 .3 GK is dominated by two low energy resonances 
just above the proton threshold (4394.9(7) keV) correspond- 
ing to two excited states in ^"S in the w 4.7 - 4.8 MeV 
range, whose 7" values were previously estimated [9] to be 3+ 
and 2+, respectively. We have observed these excited states 
in at 4688.1(4) keV and 4809.8(6) keV, respectively, via 
two separate experiments: the ^-^SipJ)^^S two-nucleon trans- 
fer reaction and an in-beam y-ray spectroscopy experiment via 
the 2^Si(3He, njf^S reaction. 

Both of our experiments result in measured resonance en- 



ergies, corresponding to the aforementioned excited states, 
which are in excellent agreement with each other Moreover, 
we have been able to reduce the uncertainty in the energies 
of these resonances with respect to what was previously ob- 
served [13] for the 3+ resonance and predicted [9] for the 2+ 
resonance. Furthermore, we have confirmed the spin-parity 
assignments of both of these resonances. As a result, our new 
^^P(;5, yj-^^S reaction rate is increased by a factor of 2 over the 
temperature range of 0.1 < J < 1.3 GK. Also, the uncertainty 
in our new rate in this temperature range has been reduced by 
72% relative to that previously determined [12]. 

This updated rate have been used to compute a full hydro- 
dynamic nova simulation which is more realistic than the post- 
processing nucleosynthesis simulation performed in Ref. [1]. 

Our new ^^P(p,y)^^S rate has only marginally improved 
the agreement between the abundances observed in presolar 
grains of potential nova origin and those obtained from sim- 
ulations. Although, our updated rate does not affect the sili- 
con isotopic abundance ratios significantly, due to a reduction 
in its uncertainty, the present nova hydro simulations can be 
compared with more reliability to the isotopic ratios measured 
in presolar grains of potential nova paternity. 

As for the nuclear structure of 3*^S, improvements in spin- 
parity assignments may be made by theoretical estimates via 
the IMME for those ^^S states whose spin-parity assignments 
are still tentative. However, this method is currently unreliable 
for A = 30 because many of the relevant analog states in ^''P 
also have unknown or tentative spin-parity assignments [40]. 
Thus, if such properties of the levels of ^^P are constrained 
better in the future, this in turn will help with the determina- 
tion of those of 3*^S. 
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